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Abstract. We construct efficient approximations for the cigcnfunctions of non-self- 
adjoint Schrodinger operators in one dimension. The same ideas also apply to the study of 
resonances of self-adjoint Schrodinger operators which have dilation analytic potentials. 
In spite of the fact that such eigenfunctions can have surprisingly complicated structures 
with multiple local maxima, we show that a suitable adaptation of the JWKB method 
is able to provide accurate global approximations to them. 
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1. Introduction 

The motivation for this study was, to a great extent, our recent paper [0, where we dis- 
covered that extreme spectral instabihty was typical for many non-self-adjoint Schrodinger 
operators, and also for complex resonances of self-adjoint Schrodinger operators with dila- 
tion analytic potentials. The full analysis of such operators involves the determination not 
only of their spectra, but also of their eigenfunctions, which may be highly non-orthogonal. 
It is especially important to have a powerful tool for approximating eigenfunctions in 
situations in which numerical computations turn out to be intrinsically unstable and, 
therefore, standard routines are insufficiently reliable. This is the case for the spectral 
problems we have been dealing with in the mentioned paper. We start by testing the 
JWKB analysis of section 2 on a typical Schrodinger operator whose spectrum has been 
studied in Then the same technique is applied to another operator with a dilation 
analytic potential. Although the two operators have different dilation analytic potentials, 
they exhibit surprisingly similar spectral properties. However, in the course of this work 
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we were primarily interested in constructing their eigenfunctions rather than comparing 
the corresponding eigenvalues. In section 3 we bring together the results of computations 
and report the conclusions drawn from their comparative analysis. This allows us to sug- 
gest that the technique we have developed can be used as a general method for a wide 
range of similar problems. 

JWKB analysis is known to be an extremely useful technique when applied to a range 
of spectral problems. To be more precise, this notion includes a vast variety of methods, 
each associated with certain asymptotic formulae. The general asymptotic theory has been 
developed in a number of books among which monographs p, |^, [1^ can be mentioned. 
There exists an extensive literature on numerical analysis of ODEs where purely numerical 
procedures are combined with and strengthened by proper JWKB methods. Examples 
of different asymptotic formulae can, for instance, be found in papers P, P, Apart 
from this, one can take advantage of similar techniques for strictly analytical purposes, 
including spectral analysis of Schrodinger operators. Recent papers on the subject include 
0, ^, 1^, 0. In the papers 0, ^ relevant JWKB-type formulae have been used to construct 
the semi-classical modes for non- self- adjoint Schrodinger operators in order to study their 
spectral properties. A major aspect of ^ is that the JWKB functions are defined 
globally, not just asymptotically at infinity, and do not involve analytic continuation to 
complex phase space. We develop the basic approach of @] in the next section, working 
out efficient eigenfunction approximations. 

Another feature of our method which makes it different from the standard approach 
was suggested by our numerical results. Having studied the spectra of several opera- 
tors, we became interested in their eigenfunctions. The latter, although more difficult 
to compute than the spectra, were also calculated by the transfer method followed by a 
conventional routine. The modes we have found are quite unexpected: some of them look 
like linear combinations of several complex Gaussian functions whose centres are quite 
distinct. It is this circumstance that prompted our decision to try to approximate the 
eigenfunctions globally by linear combinations of several JWKB approximate eigenfunc- 
tions. This permitted us to approximate rather accurately all the eigenfunctions except 
those associated with eigenvalues which are close to the origin, whereas using just one 
JWKB function usually only works for substantially larger eigenvalues. The results of 
section 3 show that in our examples each eigenvalue (starting from a certain number) 
is associated with two JWKB functions, one of them playing the leading role for lower 
and the other for higher eigenvalues. Linear combinations of these functions are proved 
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to provide fairly good approximations for the eigenfunctions involved in our examples, 
in spite of the fact that we only retain the lowest order terms in the JWKB expansion 
in this paper. Having observed this phenomenon numerically for different examples, we 
propose a method suitable for approximating eigenfunctions by several easily obtained 
JWKB functions. Their number can be arbitrary and depends on the potential of the 
operator under consideration. 

2. JWKB Formulae 

We are interested in computing the eigenvalues and eigenfunctions of the Schrodinger 
operator 

H:=-f^ + V{x) (1) 

where V{x) is a complex-valued continuous potential. Since we will adapt the JWKB 
method we study instead the Schrodinger operator acting in C2{^) given by 

dx^ 

where h > 0. Throughout the paper the potential is taken to be even; this simplification 
is not necessary, but enables us to restrict our computations to the half-line R+ instead 
of R. 

Given an eigenvalue z of Hh, let us solve the equation 

{Hf,-z)f{x) = 0. (2) 

We impose one of the conditions 

/(O) = or /'(O) = 

bearing in mind that the eigenfunctions of Hh are either even or odd. Our aim is to ap- 
proximate the true eigenfunctions f{x) by means of relevant JWKB asymptotic formulae. 
Following , we assume for now that h is sufficiently small and represent the approximate 
solution y{x) of (0) as 

y{x) = y{a + s) = ^{s; h) exp (-/i" V(-s)) (3) 

about some a > to be determined. One expects ^(s; h) to be asymptotically expanded 
in powers of h: 

oo 

i{s-h)^Y.h'Us). h^+Q. (4) 

fc=0 
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As in tp{s) is taken to be the solution of the eikonal equation 

^'^sf = V{a + s)-V{a)-r]^, ^(0) = 0, 
where real t] and a satisfy z = f]"^ + V{a). Thus, we have 

Here the sign of rj and the branch of the square root in the integrand are chosen so that 
Re ip{s) > for small s. An important consequence of the fact that V is complex- valued 
is that generically the integrand never vanishes, and thus a unique continuous branch of 
the square root is defined globally on the real line by the above formula. Substituting the 
expansion into and equating the coefficients by equal powers of h, we get a series 
of equations for C,k, k = 0,1, . . . . To normalise the function y{x) we impose the initial 
conditions {o(0) = 1, ^a:(0) = 0, /c > 0, so that y{a) = 1. A direct calculation gives us 

2eo^' + eo^"=0 

and, therefore, 

/ V{a + s)-V{a) 



-1/4 



Us) 



^2 



The function ^o(s) is globally well-defined and continuous. 

The other functions ^k, k = 1,2,..., can be successively determined as well. In this 
initial study, we omit the higher order terms in (^) and finally obtain the approximate 
solution 

yix) = eo(.)exp(-^^'eo(t)-Mt). (6) 
We substitute y directly into (0) to confirm that 

y 

Although the above analysis is only justified asymptotically as — +0, we apply the 
formulae obtained to the case h = 1, that is to the operator defined by (|I]). We refer to 
for a justification of this for large eigenvalues and for a wide class of smooth potentials. 
Thus we actually consider the eigenvalue problem 

{H-X)f = (7) 

in the rest of the paper. We expect the basic JWKB modes y{x) defined by @ to provide 
good first approximations for the eigenfunctions related to the eigenvalues A of if as 
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|A| oo. The above relation becomes the approximate formula for the eigenfunction 
of H related to A: 

y{a + s)={l j exp j dt j 

where the pair of real numbers (a, rj) solves 

X = r]^ + V{a). (9) 

Here a > since we consider x G R+; the sign of t] is chosen so that for small s in 
we have Re ip{s) > for the selected branch of the square root in the integrand. In what 
follows we shall use the key formulae (H), (H) systematically. The actual efficiency of these 
formulae can be judged by the numerical results discussed in the next section. 

Note that apart from (^) there are other ways of approximating f{x). The expression 
d^) can, for instance, be further approximated by 

This formula is easier to deal with than (P) because of its explicit expression. Naturally, 
y{x) is expected and proved to be less accurate for our aims than the original JWKB func- 
tion y{x). One can also continue the process and work out more terms ^a:, ^ = 1, 2, . . . . 
However, the formula (H) provides quite satisfactory eigenfunction approximations, which 
is justified by the results of computations related to two typical examples (see, in partic- 
ular, table 6 of subsection 3.4). 

Several final remarks are in order. The functions defined by do not generally lie in 
£2(1^); even when they are good approximations to the eigenfunctions in some interval. 
In subsection 3.2 we discuss a further truncation procedure needed to obtain £2(R') func- 
tions without sacrificing the fact that they satisfy the eigenvalue equation approximately. 
Secondly the number of solutions of (H) for a given eigenvalue A depends on the poten- 
tial V, and each of the corresponding JWKB functions may make a contribution to the 
eigenfunction, as we show in subsection 3.3. In our examples the number of solutions is 
always 0, 1 or 2. Finally the absence of turning points typical for the complex equation 
(0) should be once again mentioned. It allows us to define the functions and ^ globally 
on R without applying any special complex analysis technique. This refiects the different 
nature of the approach presented here and the standard JWKB analysis. 
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3. Numerical Experiments 



3.1. JWKB Parameters. To find the approximate eigenfunction y{x) we have to solve 
equation (P) first. In this section we consider several operators of type (0) starting by 
the harmonic oscillator operator Ho with the complex potential Vo{x) = (cx)^, c G C. In 
particular, for the m-th eigenvalue of Ho we have the equation 



c(2m + 1) = 77^ + (caf, m = 0, 1, . . . , 



which is solved exactly. The solution (a, rj) such that 77 > 0, a > defines the JWKB 
mode y{x) decaying at infinity. 

As mentioned in |I| , the eigenvalues of this and similar operators are extremely unstable 
under small perturbations, which makes numerical analysis quite difficult. To be able 
to compute higher eigenvalues and relevant eigenfunctions numerically one needs some 
additional information. For instance, when solving eigenvalue problem (|^) by the transfer 
method it is helpful to have at least rough estimates for the location of argmax|/(x)|. 
Here we take advantage of the JWKB analysis to find the centres of the eigenfunctions. 

For comparison purposes along with equation (^) let us introduce its more sophisticated 
version: 



\ + '^^{a,r^)-v'-Via) = 0. (10) 



It is obtained if we replace the eigenfunction associated with an eigenvalue A by the 
approximation 



t[x) = exp \ -ir]s H — — 



For the harmonic oscillator problem we denote the roots of by (a', r/'), those of (|T^) by 
(a", 7]"). Below we tabulate the values of r' := argmax|/(a;)| — a' and r" := argmax|/(x)| — 
a". The eigenfunctions are computed with the use of the same basic method as in |l|]. In 
our numerical exercises we take a typical value of c = e*'^/® throughout the paper. 
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Table 1. Parameters a, r for if, 



m 


a' 


a" 


argmax|/„(x) 


t' 


t" 


10 


3.371 


3.535 


3.678 


0.307 


0.143 


20 


4.711 


4.822 


4.831 


0.120 


0.009 


30 


5.746 


5.835 


5.839 


0.093 


0.004 


40 


6.621 


6.698 


6.700 


0.079 


0.002 


50 


7.393 


7.462 


7.464 


0.071 


0.002 


60 


8.092 


8.155 


8.156 


0.064 


0.001 


70 


8.736 


8.793 


8.794 


0.058 


0.001 


80 


9.335 


9.389 


9.389 


0.054 


0.0 


90 


9.897 


9.948 


9.949 


0.052 


0.001 


100 


10.430 


10.478 


10.479 


0.049 


0.001 



Equation (|I0|) easily solved by a standard iterative method provides somewhat more 
accurate information about the centres of the eigenfunctions of Ho (and other considered 
operators) than does. However, for our primary aim of eigenfunction approximation 
this is not vital, and the parameters obtained from (P) turn out to be more suitable. For 
the examples considered below we use a standard NAG routine to find the real roots of 
the equations of type and (pUj). 

An interesting example where we apply the JWKB method is an operator with a dilation 
analytic potential studied in This operator is given by 

= -c-i-^ + Vi^/^x), V{x) = exp(-a;V62)^ 
dx^ 

where cGC,6GRasin|Jl|]. Equation (|]) then becomes 

A^c = 7^ + cV{^ca) (11) 

where Am denotes the m-th eigenvalue. Another typical quantum mechanical operator 
to be considered here has a similar form to He- 

V{x) = a (exp(-7(x - Pf) + exp(-7(x + Pf)^ , a, /?, 7 G R. 

Both He and He are the operator families parametrised by a complex c; their eigenvalues 
are known to be independent of c in a sense explained in [|l|. If one puts c = 1 then the 
corresponding operators Hi and Hi are self-adjoint and known to have complex resonances 
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which are the eigenvalues of the original non-self-adjoint operators (see, for example, 
where the theory of dilation analytic resonances is exposed). 

The spectra of Ho and for sufficiently large h are proved to be very close to each 
other; the same can be said about the eigenfunctions of the two operators. For h = 100 let 
us tabulate the JWKB parameters a" found from (|TDp (where obvious changes are made 
to suit He) to make sure they also approximate the actual values of argmax|/m(a;)|. 

Table 2. Parameters a", t" for He-, b = 100 



m 


a" 


argmax|/„(a;)| 


t" 


10 


3.537 


3.682 


0.145 


20 


4.826 


4.836 


0.010 


30 


5.842 


5.850 


0.008 


40 


6.708 


6.713 


0.005 


50 


7.475 


7.479 


0.004 


60 


8.172 


8.175 


0.003 


70 


8.815 


8.817 


0.002 


80 


9.414 


9.415 


0.001 


90 


9.982 


9.983 


0.001 


100 


10.519 


10.520 


0.001 



Comparing tables 1 and 2 one can see that the results obtained for Ho and He {b = 100) 
are quite similar. The plots of the computed eigenfunctions indicate the same fact. The 
comparison of the actual eigenfunctions of these two operators contributes to our study 
of their spectra. From the numerical point of view it may turn out to be useful to have 
the values of the parameter a" found from ([To|) which approximate the centres of the 
eigenfunctions reasonably accurately. 




-30 -ao -10 10 20 30 



Figure 1. Potentials V{x) (solid line) and V{x) (dotted line); values of ^ — 
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From now on we shall concentrate on the two operators aX b = 10 and He where 
7 = 0.03, a = 100/e, /5 = 10. The choice of the parameters a, /5, 7 involved into He 
is systematic rather than random. Our idea is to take two close self-adjoint operators 
Hi and Hi and transform them into non-self-adjoint operators by the dilation analyticity 
technique (as has been done for Hi in [^). In fact, the potentials V{x) and V{x) at the 
chosen values of b and a, (3, 7 are relatively close to one another: — V"(x)||/||y(a;)|| ~ 

0.18; see also figure 1 where V, V and V — V are plotted. (Here and below by || • || 
we understand the £2-norm.) Of course, the complex potentials V{y/cx) and V{y/cx) 
thus obtained are no longer similar to each other. Remarkably, the operators He and 
He, although different, prove to have similar spectral properties. Namely, their lower 
eigenvalues are real (up to a chosen accuracy); starting from a certain number they turn 
sharply into the lower half-plane. This behaviour of the complex resonances observed in 
|I^ seems to be typical for a range of operators, and the results quoted here provide 
yet another numerical evidence of the fact. 

Below we quote the even eigenvalues of He calculated by means of the same procedure 
as proposed and implemented in our earlier paper. The eigenvalues of He are tabulated 
in [Q (see table 10 there and also figure 2 below for comparisons). 



Table 3. Eigenvalues of H, 



m 






m 






m 







4.4063 




18 


30.2243 




36 


46.9295 - 7.32702 


2 


7.3902 




20 


32.6595 




38 


48.7860 - 9.0639i 


4 


10.3846 




22 


34.8766 - 0.0005i 




40 


50.6604 - 10.9142i 


6 


13.3718 




24 


36.6986 - 0.0736i 




42 


52.5485 - 12.8733i 


8 


16.3358 




26 


38.1852 -0.6680i 




44 


54.4504 - 14.9375i 


10 


19.2609 




28 


39.8085 - 1.6721i 




46 


56.3613 - 17.1037Z 


12 


22.1312 




30 


41.5232 - 2.8670i 




48 


58.2803 - 19.3695Z 


14 


24.9295 




32 


43.2907 -4.2188i 




50 


60.2058 - 21.7326i 


16 


27.6357 




34 


45.0959 - 5.7093i 




52 


62.1367- 24.1910Z 



Given the eigenvalues, we find the values of a and t] from (pTj). In a generic situation 
this equation may have several real solutions (a^, 77^) such that > and r] has the 
appropriate sign as discussed in section 2. This is the case for both of the operators He 
and Hei unlike the former two operators. Namely, for higher eigenvalues there exist two 
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solutions providing \y{ak + s)| < \y{ak) \ for sufficiently small s. The values of a^, rj^, k = 
1,2, a2 > ai, are tabulated below (see table 4). The results for and He are qualitatively 
similar; the second solution of (0) appears for the 28-th and 30-th eigenvalue respectively. 



Table 4. Solutions of (]TT|): a, r] relate to He, a, r] to He 



m 


a 


V 


a 


fj 




m 


a 


V 


a 


V 


20 


5.0880 


3.9262 


5.4144 


3.8701 




40 


3.5726 


6.3897 


4.4905 


6.2876 


26 


5.4023 


4.3319 


5.8455 


4.2843 






12.3112 


-3.5569 


11.4371 


-3.3463 


28 


5.2223 


4.5915 


5.7562 


4.5455 




42 


3.2224 


6.6574 


4.1567 


6.5730 




10.7554 


-0.3678 










-I -I 

12.6198 


-3.9050 


11.6097 


-3.7166 


30 


5.0489 


4.8785 


5.6270 


4.8273 




44 


2.7926 


6.9310 


3.7719 


6.8549 




10.9252 


-1.3639 


10.7532 


-0.1320 






12.9812 


-4.2601 


11.7924 


-4.0723 


32 


4.7845 


5.1944 


5.4631 


5.1183 




46 


2.2971 


7.1949 


3.3172 


7.1332 




11.1804 


-1.9060 


10.8609 


-1.4247 






13.3664 


-4.6127 


11.9855 


-4.4182 


34 


4.5351 


5.5016 


5.2676 


5.4124 




48 


1.6747 


7.4508 


2.7555 


7.4085 




11.4189 


-2.3884 


10.9854 


-2.0364 






13.7846 


-4.9684 


12.1898 


-4.7583 


36 


4.2522 


5.8023 


5.0413 


5.7065 




50 


0.6079 


7.7002 


1.9911 


7.6815 




11.6864 


-2.8042 


11.1241 


-2.5247 






14.2496 


-5.3334 


12.4065 


-5.0958 


38 


3.9280 


6.0983 


4.7833 


5.9985 




52 






0.2650 


7.9532 




11.9881 


-3.1856 


11.2750 


-2.9532 










12.6375 


-5.4334 



Figure 2 explains why the second solution (a, t]) appears. Here the complex potentials 
cV{^/cx) and cV{^/cx) are plotted for < x < X, X being sufficiently large. We treat 
a; as a parameter and represent the data in the form (Re {cV), Im {cV)) for each x. The 
dots in the figure denote the values of A^c, i.e., the eigenvalues of the operators He and 
He multiplied by c (only even eigenvalues are pictured for both operators). These pictures 
refer to equation (|lTD and show how it is solved graphically. Clearly, if a dot A^c is 'inside' 
the graph then the two possible values of r]"^ = Re (A^c) — Re (cV) have opposite signs, 
while dots lying 'outside' determine two positive values of rj'^. Hence in the former case 
there is only one real pair (a,//) solving (|Tl|) (provided the proper sign of rj is chosen), 
whereas in the latter (pA]) has two solutions of this kind. As is seen, the critical points 
are m = 28 for the ffist operator and m = 30 for the second. It is from these numbers 
onwards that the behaviour of eigenf unctions changes (see figure 3 where eigenfunctions 
of He are plotted). 
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Figure 2. Potential and eigenvalues of He (left) and He (right) 

The graphs of the true eigenfunctions also suggest that for a certain range of m (28 < 
m < 40 and 30 < m < 44 for H^. and H^. respectively) the solution of (|^) has to be sought 
as a linear combination of two JWKB modes rather than a single function of type (^. 
These modes are determined by the corresponding values of (a, rj). Of course the problem 
d^) itself can have only one £2 eigenfunction for each eigenvalue, but the peculiarities of 
the solutions indicate that the existence of more than one JWKB approximation has great 
significance for our problem. 

3.2. Cut-Off Technique. There is a possibility that Re {—i^S,o{s)~'^) > for some 
range of s for particular JWKB modes. Therefore, we cannot use the formula (P) to 
approximate £2-functions / on the whole R+. In fact, most of the modes defined by 
d^) in our examples, although good around their centres where s is small enough, either 
increase as s — > cxd or change rapidly near s = —a. 

When this problem occurs, we modify the original functions y defined by (|^) putting 

yi[x) yi[si) y2[x) y2{s2) 

where the points si, S2 are chosen properly for each JWKB function. Otherwise the 
functions are not changed anywhere. Our aim is to replace the formerly defined modes by 
decaying exponentials beyond [s2,Si] if necessary; we shall keep the same notation y{x) 
for the modified functions in what follows. 

It is clear that our results depend on the choice of these Sk- Generally, in practical 
computations we have to vary two parameters to minimise the approximation error for 
each particular eigenfunction (see the next subsection for corresponding formulae). To 
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get optimal results one can take the initial values of Si and S2 reasonably close to the 
corresponding local minima of \yi\ and I1/2I respectively, replacing by exponentially 
decaying functions according to ([l2|) . Then the values of Sk providing the best possi- 
ble approximations among those given by the above formulae are found by a standard 
minimisation procedure. 

However, this general routine can be substituted by a simple although not quite univer- 
sal method. We observe that in the examples considered here a possible way of choosing 
Sk is as follows. For the original functions yk we evaluate 



Sk = argmax 



Re 



then take 

Vkix) 

The reasoning leading to this choice, although not rigorous, may prove to work in a generic 
situation. Replacing y'k/yk by constants beyond the points Sk as in (JT^), we obtain our final 
modes which are almost orthogonal to one another. Their suitable linear combinations 
approximate the actual eigenf unctions of our operators fairly well. This is the subject of 
the remaining subsections. 

3.3. Constructing Eigenfunctions. In a generic situation when equation (|^) has sev- 
eral real roots of proper signs, we have to approximate the relevant eigenfunction by the 
linear combination of the corresponding JWKB functions: 

n 

(l)ix) = Y'^kykix). (13) 

k=l 

Here yk{x) are basic JWKB modes related to (a^, rjk). The number of modes entering ( [TB| ) 
can be arbitrary; in our examples n = 1 or 2. The coefficients Ck are to be determined. 
One approach to computing these constants is as follows. Let minimise the value of 

^ ^ ||/(x)-0(x)|| ^ ^^^^ 



Denote d{x) := f\x) — 4>{x)] a direct calculation gives us 



11/(3;) - Ckyk{x)f = Wff - 2Re Y ^kUk + Y (^jkCkCj I (15) 

k=l \k=l 




where 



00 



'^k = J fVkdx, ttjk = J ykVjdx. 
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Differentiating (plSf ) with respect to and putting the variation equal to zero, we obtain 
the desired constants: 



C = A ^U, C = {ci, . . . ,CnY , U = {ui,. . . ,Un)^, A={ajk),j,k 



,n. 



(16) 



The integrals involved into ([Tq ) are easily computed by means of the routine described in 
|jl|] in terms of /'//, y't/Uk- This procedure does not require computing the eigenf unctions 



themselves: it is only the values of m^, a,fc, j, k 



n, that we have to calculate. Then 



the constants c^, k = 1, . . . ,n, are found from (16). Obviously, for a given eigenfunction 
our Ci, C2 must depend on the normalisation of the functions / and yk while the value of 
|c2/ci| must not. We take ||/|| = ||ei|| = ||e2|| = 1 to make the matrix A well-conditioned 
and then repeat computations for several different values of ||efc|| to double-check our 
results (we leave ||/|| = 1 throughout the paper). The procedure developed in |1| turns 
out to be suitable for our purposes and provides reliable answers. 

Let us tabulate c := |c2/ci| for the examined eigenvalues of the operators He and He 
(as above, m denotes the number of an eigenvalue). 

Table 5. Values of c for the operators He and He 



m 


He 


He 




m 


He 


He 


28 


0.0447 






42 


277.2584 


13.5057 


30 


0.0928 


0.0328 




44 


1434.62 


48.2824 


32 


0.2083 


0.0621 




46 


7615.34 


177.1347 


34 


0.7876 


0.1420 




48 


41623.37 


660.8588 


36 


2.9866 


0.3890 




50 


354757.06 


2486.87 


38 


13.4820 


1.1897 




52 




10590.02 


40 


65.6680 


3.9086 











The contents of the above table is illustrated qualitatively by the plots of the eigenfunc- 
tions (see figure 3). Indeed, for a certain range of eigenvalues the relevant eigenfunctions 
have two distinct components clearly seen in the plots. For lower eigenvalues the first 
mode yi is dominating, while the contribution of y2 is more significant for the high energy 
spectrum. All our numerical results including the plots of the eigenfunctions indicate the 
transition from the first leading JWKB function yi to ?/2 in formula (|13]). This phenome- 
non is observed for both examples studied here. 
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In the previous subsection we mentioned the parameters entering the modified func- 
tions Hk- The formulae of this subsection contain these parameters fixed at some starting 
values. After minimising A with respect to Ck as described we minimise the obtained 
function with respect to Sk, and thus get the final results to be tabulated in the next 
subsection. 



3.4. The Efficiency of Approximations. Having found c^, = 1, 2, we then compute 
A defined by (|1^ in terms of the calculated coefficients Uk, ajk- The values of A tabulated 
below characterise the accuracy of the JWKB approximations we use. 

As is seen from table 6, not only do the JWKB modes approximate eigenfunctions 
corresponding to higher eigenvalues but they also turn out to be surprisingly efficient 
even for the lower part of the spectrum. The plots of \dm\ in comparison with \fm\ (see 
figure 4) also indicate the effectiveness of the approximations obtained by our approach. 
The corrections dm for other values of m either look similar to those shown in figure 4 or 
are invisible compared to the corresponding eigenfunctions. 

These numerical results complete our study of eigenfunction approximations. 



Table 6. Minimal values of A 



m 


He 


He 




m 


He 


He 


10 


0.132627 


0.153114 




32 


0.015264 


0.015016 


12 


0.042214 


0.048187 




34 


0.017776 


0.015332 


14 


0.029275 


0.032804 




36 


0.014177 


0.020596 


16 


0.021119 


0.024966 




38 


0.009434 


0.025885 


18 


0.015969 


0.019302 




40 


0.006557 


0.025631 


20 


0.012366 


0.011576 




42 


0.004796 


0.012372 


22 


0.011045 


0.009765 




44 


0.003873 


0.007280 


24 


0.010724 


0.009110 




46 


0.002828 


0.005568 


26 


0.014866 


0.008487 




48 


0.002236 


0.004359 


28 


0.015133 


0.009798 




50 


0.001732 


0.003464 


30 


0.015232 


0.012214 




52 




0.003133 
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2 4 6 6 10 12 U 2 4 6 8 10 12 14 




5 10 15 20 5 10 15 20 




Figure 3. Eigenfunctions of He- Re fm, ^ai f^, \fm\- Prom left to right, top to bottom: 

m = 28, 30, 32, 34, 36, 38 
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